Shock Acceleration and 7-Ray Emitting Supernova Remnants 



Matthew G. Baring 1,2 , Donald C. Ellison 3 , Stephen P. Reynolds 3 , 
Isabelle A. Grenier 4 and Philippe Goret 4 

1 LHEA, NASA/Goddard Space Flight Center, Greenbelt, MD 20770, USA 
2 Compton Fellow, Universities Space Research Association 
^Department of Physics, North Carolina State University, Raleigh, NC 27695, USA 
4 DSM/DAPNIA / Service dAstrophysique, CE-Saclay, 91191 Gif-sur-Yvette, France 

Abstract 

Diffusive shock acceleration in the environs of a remnant's expanding shell is a popu- 
lar candidate for the origin of SNR gamma-rays, as well as providing the principal source 
of galactic cosmic rays. In this paper, results from our study of non- linear effects in 
shock acceleration theory and their impact on the gamma-ray spectra of SNRs are pre- 
sented. These effects describe the dynamical influence of the accelerated cosmic rays on 
the shocked plasma at the same time as addressing how the non-uniformities in the fluid 
flow force the distribution of the cosmic rays to deviate from pure power-laws. Such devi- 
ations are crucial to gamma-ray spectral determination. Our self-consistent Monte Carlo 
approach to shock acceleration is used to predict ion and electron distributions that spawn 
neutral pion decay, bremsstrahlung and inverse Compton emission components for SNRs. 
We demonstrate how the spatial and temporal limitations imposed by the expanding SNR 
shell quench acceleration above critical energies in the 500 GeV - 10 TeV range, thereby 
spawning gamma-ray spectral cutoffs that are quite consistent with Whipple's TeV upper 
limits to the EGRET unidentified sources that have SNR associations. We also discuss the 
role of electron injection in shocks and its impact on the significance of electromagnetic 
components to GeV-TeV spectral formation. 



1. Introduction 

Supernova remnants have long been invoked as a principal source of galactic cosmic rays, cre- 
ated via the process of diffusive Fermi acceleration at their expanding shock fronts (e.g. Drury 
1983, Lagage and Cesarsky 1983). They can also provide gamma-ray emission via the inter- 
action of the cosmic ray population with the remnant environment; this concept was explored 
recently by Drury, Aharonian and Volk (1994). In their model, the gamma-ray luminosity is 
spawned by collisions between the cosmic rays and nuclei from the ambient SNR environment. 
In the more recent models of Mastichiadis and de Jager (1996), Gaisser, Protheroe and Stanev 
(1997), and Sturner, et al. (1997), ee and ep bremsstrahlung, and inverse Compton scattering 
involving shock-accelerated electrons interacting with the cosmic microwave background and 
also IR/optical emission (from dust/starlight) form added components. See the reviews by de 
Jager and Baring (1997) and Volk (1997, this volume) for a discussion of various models. 

With no definitive detections of gamma-rays from known supernova remnants, the moti- 
vation for modelling these "hypothetical" sources hinges on a handful of spatial associations of 
unidentified EGRET sources (Esposito et al. 1996) at moderately low galactic latitudes with 



well-studied radio and X-ray SNRs. These include IC 443, 7 Cygni and W44, and most have 
neighbouring dense environments that seem necessary in order to provide sufficient gamma-ray 
luminosity to exceed EGRET's sensitivity threshold. The remnants associated with several of 
the EGRET unidentified sources show an apparent absence of TeV emission, as determined by 
Whipple (Lessard et al. 1995, Buckley et al. 1997), which could be explained by an intrinsic 
cutoff in the SNR-generated cosmic ray distribution (e.g. Mastichiadis and de Jager 1996). 

All of the above models invoke simple power-law accelerated particle populations. In 
this paper, we utilize the more sophisticated output of shock acceleration simulations (e.g. 
Jones and Ellison 1991) to address the issues of spectral curvature and the maximum energy 
of acceleration in the context of SNR gamma-ray emission. We use output from the fully non- 
linear, steady-state Monte Carlo simulations of Ellison, Baring and Jones (1996) to describe 
the accelerated particles in environments where they influence the dynamics of the SNR shell. 
Our results make clear predictions of what maximum energies of gamma-rays are expected and 
more accurate predictions of the level of TeV emission in these sources. Our model provides a 
prescription for defining realistic values of the non-thermal electron/proton abundance ratio. 

2. Fermi Acceleration at SNR Shocks 

Shock acceleration is usually assumed in astrophysical models to generate power-law particle 
populations. This approximation omits the effect the accelerated particles themselves have on 
the hydrodynamics of their shocked environment. Such non-linear effects, well-documented in 
the reviews of Drury (1983) and Jones and Ellison (1991), have a feedback on the acceleration 
mechanism and its efficiency. The slope of the cosmic ray distributions depends purely on 
the compression ratio r = Wups/wdown °f fl° w speeds either side of the shock, and this ratio 
ultimately depends on the shape of the ion distributions. Our kinematic Monte Carlo technique 
(see Ellison, Baring and Jones 1996) for simulation of Fermi acceleration in the non-linear 
regime is ideal for fully exploring particle populations, acceleration efficiencies and spectral 
properties. The simulation follows particle convection and diffusion in the shock environs using 
a simple scattering law: a particle's mean free path A is proportional to its gyroradius r g . 
Generalizations of this law are easily handled by the Monte Carlo simulation. 

Clearly the maximum energy of Fermi-accelerated ions is of central importance to TeV 
observations. It can be determined in the Sedov phase by equating the acceleration time to the 
remnant age t SNR , producing values of the order of (see Baring et al. 1997): 

5- ~ 4.7 Z=l 2 (Bj,) (ig^) 2 ^ TeV (1) 

where B x is the field in units of Gauss, V s ^ is the expansion speed, and rj — X/r g . In the Sedov 
phase, where the peak luminosity is expected (see Drury, Aharonian and Volk 1994), E MAX is 
actually a slowly increasing function of time (see Baring et al. 1997). For young SNRs early 
in the Sedov phase ( 200- 5000 yrs), typically ~ 200 -3000 km/sec, implying age-limited 
acceleration terminating at around 1 — 10 TeV. The diffusion scale of the SNR medium is k/u 
for a diffusion coefficient of k — Xv/3 , so that Eq. (1) generally guarantees that maximum 
diffusion lengths are always considerably smaller than the shock radius. 

Typical simulation output is shown in Fig. 1, where the resulting ion and e~ distributions 
are exhibited. The pressure of the accelerated ions (p and He 2+ ) acts to slow down the fast- 
moving flow upstream of the shock, creating a maximum compression ratio r that is much 
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Fig. 1: Proton, He 2+ and e~ distributions resulting from the Monte Carlo simulation of accel- 
eration at SNR shocks at three different ages in the Sedov phase. The transition to relativistic 
energies produces "bumps" in the spectra (e.g. at around 1 GeV for protons). The ion spec- 
tra are the upper two histograms at GeV-TeV energies (protons - solid, helium - dashed), 
providing nearly all of the total energy density of the system. The dashed-dot histograms rep- 
resent the e~ distributions (see next section for a discussion of the parameters E CT - lt and the 
injection temperature T e in j ). As the shock slows with age, it weakens (r declines). Upstream 
parameters were a density of ni — 1 cm -3 and field strength B\ = 3/x G. 

greater than the canonical "strong shock" value of 4. Such large r are realized at the largest 
scales, implying an increased "efficiency" of accelerating higher energy particles than those at 
lower energies. Hence upward curvature appears (e.g. see Baring et al. 1997) in the proton, 
a and electron distributions: high energy particles have longer mean-free paths and therefore 
typically influence the flow on larger scalelengths. The escape of particles upstream of the 
shock (i.e. outside of the remnant's shell), which occurs on a timescale of ~ t SNK , gives a 
cessation of acceleration, in this case at TeV energies, forcing the compression ratio to increase 
to compensate in the flow hydrodynamics. These non-linear predictions of the spectrum differ 
significantly from the test-particle case (Ellison, Baring and Jones 1996), with deviations from 
power-law behaviour by as much as a factor of around 3 over four decades in particle energy: 
this is a significant influence on predictions of TeV fluxes in gamma-ray SNRs. 

3. Gamma-Ray Emission from Shell- Type Remnants 

In applying the full non-linear Monte Carlo simulation to the prediction of SNR 7-ray emission, 
we model the tt° decay component, bremsstrahlung and inverse Compton scattering of back- 
ground radiation fields. Cosmic ray ions collide with nuclei in the cold ambient ISM to produce 
7T° s, which subsequently decay to create two photons; the decay spectra were calculated much 
along the lines of the work of Dermer (1986). Bremsstrahlung and inverse Compton contribu- 
tions were calculated along the standard lines used in Gaisser, Protheroe and Stanev (1998) 



and Sturner et al. (1997). Fig. 2 shows the emergent 7-ray spectra resulting from particle 
distributions (similar to those in Fig. 1) generated by the non-linear Monte Carlo simulation. 
The non-linear modifications can alter the TeV/EGRET flux ratio by as much as a factor of 
2-3 (Baring, Ellison and Grenier 1997). In this case, ions are accelerated out to around a few 
TeV per nucleon, providing compatibility with the Whipple upper limits. Turnovers at such 
low energies are not predicted in the work of Gaisser, Protheroe and Stanev (1998), but are in- 
voked by Mastichiadis and de Jager (1996), and Sturner et al. (1997). One implication of such 
low maximum energies for cosmic rays may be that gamma-ray SNRs are a gamma-ray bright 
minority of the remnant population, with other SNRs being required to produce the bulk of 
cosmic rays out to the 10 14 eV "knee." Validity of such a contention would be contingent upon 
a confirmation that at least some of the EGRET unidentified sources in Esposito et al. (1996) 
represent detections of shell-related emission. Alternatively, such emission may be common 
but at flux levels below EGRET's sensitivity, awaiting future detections; such a scenario might 
permit most shell-type 7-ray remnants to simultaneously be prolific producers of cosmic rays. 




-15 -12 -9-6-3 3 6 
Log 10 E r [MeV] 

Fig. 2: Inverse Compton (IC), bremsstrahlung and pion decay (p-p) 7-ray emission (Hux) 
spectral components, and radio to X-ray synchrotron radiation, obtained by integrating over 
shock acceleration-produced ion and electron distributions much like those in Fig. 1. These 
were computed assuming a source volume of 1 pc 3 and distance of 1 kpc, and the shock model 
had parameters u nps = 876 km/sec, n Pi i = 0.25 cm -3 and field strength B\ = 3/i G. Such a 
model has insufficient density to compare with the EGRET data for the source associated with 
IC 443 (Esposito et al. 1996) that are plotted, so we depict via the heavy solid curve a scaling 
of the total spectrum by a factor of 500; this could be realized through a larger emission volume 
or source distance. Also plotted are radio data from the compilation of Erickson and Mahoney 
(1985), the Ginga Hux (Wang et al. 1992) and upper limits from the Whipple (Buckley et al. 
1997) and HEGRA (Prosch et al. 1996) experiments. 



Some features of the model results illustrated in Fig. 2 together with a collection of data 
for IC 443 that are immediately apparent are that (i) the low density chosen (following Gaisser, 
Protheroe and Stanev 1998) inhibits pion decay emission, compared with the high densities 
inferred by Drury, Aharonian and Volk (1994) in matching EGRET flux sensitivities, (ii) the 
strong contribution of the bremsstrahlung in hard gamma-rays generates a spectrum steep 
enough to be compatible with the EGRET source spectral index, and (hi) the 1-10 TeV-range 
turnovers in the particle distributions provide total consistency of the radiation spectrum with 
Whipple and HEGRA upper limits. This last feature refutes recent suggestions that the upper 
limits obtained by Whipple to IC 443 are so constraining that they cause serious problems 
for the Fermi mechanism (see Baring 1997, and de Jager and Baring 1997, for more detailed 
discussions). Scaling up the source distance and volume in the figure emphasizes the suitability 
of this particular model for the gamma-rays, but provides inconsistencies with the radio and 
Ginga X-ray data. Hence shock acceleration can comfortably accommodate the gamma-ray data 
points, and modification of SNR parameters to reduce the relative importance of bremsstrahlung 
can yield compatibility with the X-ray data (discussed in Baring et al. 1997). The flat radio 
synchrotron spectral index poses more of a problem (see Baring et al. 1997). We note that the 
localization of EGRET gamma-rays from IC 443 in Esposito et al. (1996) is entirely consistent 
with a shell origin, in contrast to the more confining localizations available for 7-Cygni (Brazier 
et al. 1996, discussed in de Jager and Baring 1997). For this reason, we view IC 443 as a better 
candidate for comparing models with gamma-ray data. 

One important aspect of our work is the connection we make between properties of thermal 
and suprathermal electrons and the resulting e/p ratio. We inject electrons as "test particles" 
at temperature T ei i n j , which is usually in the keV range (see Fig. 1) so as (i) to be compatible 
with temperatures deduced from X-ray observations, and (ii) to mimic a shock dissipation that 
deposits a sizeable fraction of the proton Rankine-Hugoniot thermal energy into the electron 
population. Furthermore, we define a critical energy E cv - lt , above which electrons efficiently 
resonate with Alfven- whistler modes (with A oc r g ), and below which electron diffusion is 
inefficient, as is expected at suprathermal energies due to wave damping. These choices define a 
range of possibilities that describe how electrons are injected into the Fermi process. Electrons 
"injected" at lower energies (i.e. lower T ejin j and E CTit ) have shorter mean free paths, and 
consequently steeper spectra. Hence these parameters determine the value of the e/p ratio 
at 1 GeV (and above) which was around 0.7 for the case exhibited in Fig. 2 (which had 
2e,mj = 9.4 x 10 5 K and E CV i t = 100 keV). Generally, our models produce (see Baring et al. 
1997) e/p ratios considerably less than unity (e.g. see Fig. 1), which are consistent with cosmic 
ray abundances (e.g. Miiller et al. 1995) and modelling of the diffuse 7-ray background (Hunter 
et al. 1997), but provide worse fits of the gamma-ray data for IC 443 that are depicted in Fig. 2. 
The shape of the gamma-ray spectrum below 1 GeV is strongly dependent on the e/p ratio 
above 1 GeV. We believe that electron diffusion in turbulent plasmas favours lower injection 
energies. Consequently, we contend that lower levels (than in Fig. 2) of inverse Compton 
and bremsstrahlung emission relative to pion decay radiation are far more likely (contrasting 
Sturner et al. 1997) to be representative of 7-rays from SNRs. Notwithstanding, for electrons 
and protons accelerated to similar energies above 30 TeV, inverse Compton scattering is a more 
efficient creator of super- TeV gamma-rays, and is expected (Mastichiadis and de Jager 1996) 
to dominate such signals from sources like SN1006. The recent detections of SN1006 by the 
CANGAROO experiment above 1.2 TeV, announced at this meeting (see Tanimori et al., these 
proceedings), have ushered in a new era for supernova remnant studies. 



4. Conclusion 



In conclusion, the non-linear effects addressed here are desirable input to any model that 
invokes Fermi acceleration at SNRs. They produce cutoffs in the spectrum that can comfortably 
accommodate Whipple's TeV upper limits to unidentified EGRET sources, and provide a range 
of realistic estimates of the non-thermal e/p ratio. It is clear that if some of the EGRET 
detections turn out to be of gamma-rays generated in the environs of remnant shells, then 
gamma-ray emitters must be a minority of remnants, perhaps mostly young, given that they 
cannot produce ions above around a few TeV in profusion. Remnants that provide cosmic rays 
up to the knee must consequently be a gamma-ray quiet majority. Alternatively, if fluxes of shell 
origin are well below EGRET's and Whipple's flux sensitivities, then the notion that shell-type 
remnants are simultaneously gamma-ray bright and prolific producers of cosmic rays becomes 
tenable. It has therefore become evident that the Whipple and HEGRA upper limits have not 
destroyed the hypothesis that shocks in shell-type remnants energize the particles responsible 
for the gamma-ray emission, but rather have provided a powerful tool for constraining our 
understanding. Given the recent detections of SN1006 by CANGAROO, we expect that in the 
near future, coupled TeV/sub-GeV, MeV and X-ray observations will discriminate between the 
various models, and refine our understanding of gamma-ray SNRs. 
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